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ABSTRACT 

We present new light curves covering 14 to 19 years of observations of four bright proto¬ 
planetary nebulae (PPNs), ah O-rich and of F spectral type. They each display cyclical light 
curves with significant variations in amplitude. All four were previously known to vary in light. 
Our data were combined with published data and searched for periodicity. The results are 
as follows: IRAS 19475+3119 (HD 331319; 41.0 days), 17436+5003 (HD 161796; 45.2 days), 
19386+0155 (101.8 days), and 18095+2704 (113.3 days). The two longer periods are in agreement 
with previous studies while the two shorter periods each reveal for the first time reveal a dominant 
period over these long observing intervals. Multiple periods were also found for each object. The 
secondary periods were all close to the dominant periods, with P 2 /Pi ranging from 0.86 to 1.06. 
The variations in color reveal maximum variations in T e ff of 400 to 770 K. These variations are 
due to pulsations in these post-AGB objects. Maximum seasonal light variations a re all less t han 
0.23 mag (V), consistent for their temperatures and periods with the results of Hrivnak et al.1 
( 2010h for 12 C-rich PPNs. For all of these PPNs, there is an inverse relationship between period 
and temperature; however, there is a suggestion that the period—temperature relationship may 
be somewhat steeper for the O-rich than for the C-rich PPNs. 


Subject headings: planetary nebulae: general — stars: AGB and post-AGB — stars: oscillations — stars: 
variable: general 


1. INTRODUCTION 


Proto-planetary nebulae (PPNs) are objects in 
the short-lived evolutionary phase between asymp¬ 
totic giant branch (AGB) stars and planetary neb¬ 
ulae (PNs). This transitional phase in the evolu¬ 
tion of low- and intermediate-mass (0.8—8 Mq) 
stars is expected to be several thousand years 
long 1 depending upon the mass of the star ( Blocked 
19951 1. This relatively short time scale contributed 


to the initial difficulty in identifying stars in this 
stage prior to the availability of data from the In- 
fraRed Astronomical Satellite (IRAS), which be¬ 
came available in 1984. IRAS allowed the de- 


1 Present address: Adler Planetarium, 1300 

S. Lake Shore Drive, Chicago, IL 60605, USA; 
knault@adlerplanetarium.org 


tection of candidate objects based on their mid- 
infrared emission from circumstellar dust formed 
during the late-AGB phase. A good sum mary of 
their discovery and properties is given by iKwokl 
( 2000l ). while their place in the br oader context 
of post-AGB stars is described by Van Winckel 
(1200311 . 


It soon became apparent that some of the PPNs 
were oxygen-rich (O-rich) and some were carbon- 
rich (C-rich). This was initially based on their 
mid-infrared spectra as obtained by IRAS and 
later by the Infrared Space Observatory (ISO). 
Those that were 0-ricli showed amorphous sili¬ 
cate features at 9.7 and 1 8 /im, bot h in emis¬ 
sion and absorption ( Volk et ahl 199llf . and ISO 
observations also revealed crystalline silicates at 
a number of wavelengths between 20 and 43 p,m 
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(jMo lster et all 120021 ). Those that were C-rich in¬ 
stead showed infrared aromatic b ands between 3 
and 12 //in ( Hrivnak et al.l 12000 ) that are usu¬ 
ally identified as due to polycyclic aromatic hy¬ 
drocarbon (PAH) molecules. High-resolution vis¬ 
ible spectra supported these mid-infrared chem¬ 
ical classifications, with the O-rich showing C 
to O ratios (C/O) of ~0.5 (summarized in this 
pape r) and the C-rich showing C/O ratios of 
1—3 ( Van Winckel fe Revniersll2000t iReddv et al. 
20021 ). The signature of the O-rich or C-rich 


nature can also be seen in the gaseous circum- 
stellar nebulae, which can show OH (O-rich) or 
HCN (C-rich) emission in the centim eter an d mil¬ 
limeter radio regio ns, respectively dLik kcl 1989; 
Omont et al.lll993 ). 

Photometric observations of the brigh ter PPNs 
show e d that they varied in brightnes s (Arkhi pova et al 
2000l : Hrivnak &; Lu 20001 : Hrivnak et al.1 20021) . 
Subsequent long-term photometric monitoring has 
revealed period s for many of them, ranging from 
35 to 160 days dArkhipova et al.1 2000L 2010 , 2011 


Hrivnak et al. 2010h . When combined wit h stellar 
pulsation models (eg., iFokin et~ah 2001 ). pulsa¬ 
tion periods and amplitudes have the potential of 
yielding information on the mass, luminosity, and 
stellar structure of these transitional objects. 

In this paper, we analyze the light curve vari¬ 
ations of four bright O-rich PPNs. Our obser¬ 
vations range from 1994 to 2012. We also in¬ 
clude in our analysis the published complemen¬ 
tary photometric observations of these four to 
increase the data sample. This combined sam¬ 
ple extends the temporal coverage, fills in gaps 
in the data, and increases the density coverage. 
In these ways it provides an excellent sample 
from which to search for multiple periods and to 
compare with the res ults of our previous ly pub¬ 
lished C-rich sample ( Hrivnak et al . m' Paper 
I). There is observational evidence that the tran¬ 
sition from 0-ricli to C-rich chemistry occurs for 
stars with initial masses o f approximately 1.5 to 


stars witn initial masses ot approximately 1.0 to 
3 Mq ( Lattanzio &; Woo'd! 20041) . Thus one might 
expect different pulsational properties between O- 
rich and C-rich PPNs, since they result from stars 
of different masses and luminosities. 


2. PROGRAM OBJECTS 

The four program objects are listed in Table [T] 
along with some basic information. The spectra 
of all four have been classified as F supergiants. 
All have high-resolution spectral abundance stud¬ 
ies which not only show that they are O-rich, 
but that they are iron-poor, with [Fe/H] rang¬ 
ing from —0.25 to —1.1. Three of them display 
double-peaked spectral ene rgy distributions (SED; 
Hrivnak et al. 1988 . 19891) . This is characteristic 
of PPNs, with the first peak due to the light of the 
central star, dimmed and reddened by circumstel- 
lar dust, and the second peak due to the thermal 
emission from this dust. IRAS 19386+0155 dis¬ 
plays only a single peak in its SED. Below is a 
brief background description of each. 

2.1. IRAS 17436+5003 

Initial attention to the unusual properties of 
IRAS 174 36+5 003 (HD 161796, V814 Her) was 
drawn by iBidelmanl (Il98ll) . who listed it among 
a small group of four supergiants located at high 
galactic latit ude. FbHowingthe launch of IRAS , 
Partlrasarathv fc Pot tasch ( 19861) noted its large 
infrared excess and sugge sted that it m i ght b e 
a post-AGB object, as did iHrivnak et al.l ( 19891) . 
who included it in a multi-wavelength study of 
a eight PPN candidates of F—G spectral types. 
Imaging with the Hubble Space Telescope ( HST) 
revealed a small elliptical nebula with size o f 4.4" 
_x_2.5" surrounding the bright star ( Ueta et all 
2000h . The mid-infared ISO spectrum shows 
the broad amorphous silicate emission features at 
9.7 and 18 pm and additional crys talline silicate 


featu res in the 15—70 pm region (jMolste r et al 
20021 ) . 


2.2. IRAS 18095+2704 

The basic properties of IRAS 18095+2 704 
(V887 Her) were initially discussed bv lHrivnak et al 
( 19881) . who identified the IRAS source with a 
then un-cataloged 10th magnitude star. Its HST 
image shows an elliptical or perhaps bipolar neb¬ 
ula with a size of 1.9" x 1. 2" surrounding the 
bright star ( Ueta et al.l 2000l) . The mid-infared 


ISO spectrum displays silicate emission at 9.7 and 
18 pm , and weak emission features in the region 
from 10 —27 pm that are attrib uted to crystalline 
silicates ( Hrivnak et al. 2015bl) . This object is an 
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OH maser llEder et alJIlQSi OH 53.8+20.1). 


2.3. IRAS 19386+0155 

The basic properties of IRAS 19386 +0155 


(V1648 Aql) were initially discussed by van der Veen, Habing} l fe J fceWa(le^ X 
, nonK • . j r . a ^ 11 > . , ^ , witii cacii to match the stanc 


( 1989h in their study of pos t-AGB transitional ob¬ 
jects. iPereira et al.l (120041 1 carried out a visible 
and infrared spectroscopic study to determine the 
chemical composition of the star and to model its 
envelope. The SED shows only one broad peak 
reaching maximum energy in the region from 4.5 
to 25 pm. This seems to indicate a relatively ob¬ 
scured and reddened central star. It has not been 
observed with HST, but it has been observed and 
resolved in the mid-infrared, with a size of 2 . 6 " 
and possessing a br ight core and extended halo 
( Lagadec et al.ll201lh . The mid-infrared ISO spec¬ 
trum shows the 18 pm silicate emission feature and 
features in the 10 pm region that we interpret as 
self-absorption in the middle of a silicate emis¬ 
sion feature. ( Hrivnak et al. 20151^ OH maser 


emission has been detected (jLewisI 2000 ). 


2.4. IRAS 19475+3119 

The basic properties of IRAS 19475+3119 (HD 
331319, V251 3 Cyg) were broadly discussed by 


Garcia-Lario et al. (jl997 ). HST imaging shows 


that the nebula has an interesting, quadrupolar 


shape surrounded by a faint halo dSahai et al 
2007 ). The mid-infrared ISO spectrum shows 


some features identif ied as crystalline silicates 
( Sarkar fc Sahail 2006). 


3. PHOTOMETRIC OBSERVATIONS 

Photometric observations were carried out at 
the Valparaiso University Observatory (VUO) us¬ 
ing the 0.4-rn telescope and a CCD detector. They 
began in 1994 and continue through the present 
time; those reported in this study run through 
2012, except for IRAS 19386+0155, for which we 
stopped observations in 2007. Initial observations 
were made primarily through the V filter, with 
some occasionally through the Rc filter, but the 
frequency of Rc observations increased and since 

2000 they they have been made together. Unfortu¬ 
nately a problem developed with the V filter which 
caused us to reject the V data from the 2000 and 

2001 seasons and part of 2002. This led to a gap 
in each of the Plight curves. Beginning in 2007, B 


observations were also made for all of the sources 
except IRAS 19386+0155. In 2008 we upgraded 
from a Photometric Star I CCD, with a field of 
view of 816 x 517, to an SBIG 6303 CCD, with a 

Filter sets were used 
standard Johnson B and V 
and Cousins Rc systems. The observations were 
reduced with IRAF0, using standard procedures 
to remove cosmic rays effects, subtract the bias, 
and flat field the images. Aperture photometry 
was carried out with an aperture of 11 " diameter. 

Given the frequently partially-cloudy nights at 
our site, we carried out a program of differential 
photometry. Three comparison stars were used 
for each of the program stars. Ideal comparison 
stars would be a good match in color and bright¬ 
ness, but that was not achievable in these cases 
due to the brightness of the program stars and 
the size of the field of view of our CCD detec¬ 
tors. The comparison stars used are listed in Ta- 
bleEl When we changed to the larger format CCD, 
we were able to include new, brighter comparison 
stars C 2 and C 3 for IRAS 17436+5003, and these 
are the ones that we have listed in Table El For 
each PPN, the constancy of the primary compar¬ 
ison star, Ci, was examined with respect to the 
other two comparison stars. The precision of the 
observations is not as high as we would like due 
to the large differences in brightness between the 
program and comparison stars. This restricted the 
integration times to avoid non-linearity in the pro¬ 
gram stars, but consequently lead to less precision 
in the comparison star observations. We find the 
comparison stars to be reasonably constant over 
time, to within ±0.015 mag, much smaller than 
the variations found in the program stars. Com¬ 
ments about individual comparison stars are found 
in the footnotes to Table El In Table [3] are listed 
the numbers of observations of the program stars 
in each filter, along with the statistical uncertain¬ 
ties in the differential magnitudes. The typical 
precision in the differential magnitudes is ±0.006 
mag. 

Standardized photometry of these program star 
fields was carried out on two different nights at the 
VUO, and the standard magnitudes are listed in 


1 IRAF is distributed by the National Optical Astronomical 
Observatory, operated by the Association for Universities 
for Research in Astronomy, Inc., under contract with the 
National Science Foundation. 
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Table [2] The differential magnitudes of the PPNs 
were also transformed from the instrumental to 
the standard magnitude systems. Standardized V 
and Rq differential photometry for the four pro¬ 
gram stars from 1994 to 2007 is listed in Table [IJ 
except for the 2007 observing sets that included B 
also. Those are listed in Table [5] along with BVRc 
observations made with the newer CCD beginning 
in 2008. 

4. VUO LIGHT CURVES 
4.1. IRAS 17436+5003 

The VRc light and color curves of IRAS 
17436+5003 are displayed in Figure |T] The light 
curves possess a total peak-to-peak variation of 
0.23 mag in V and 0.18 mag in Rc, with a rel¬ 
atively large variation in the yearly amplitudes 
and a slight variation in the yearly means. This 
variability is superimposed perhaps on a slight, 
non-monotonic trend of increasing brightness of 
~0.04 mag (V) over the 19 years of our obser¬ 
vations. The peak-to-peak range within a season 
varies greatly, from 0.05 (2004) to 0.21 (2010) mag 
in V and from 0.06 (2004) to 0.18 (2010) mag in 
Rc, with an average ratio of the seasonal range of 
Rc to V of 0.90. The variations are even larger 
in B and (B— V) , with an average ratio of the sea¬ 
sonal B to V of 1.36 over the six years in which we 
have B observations. Samples of the B light and 
the ( B— V) color curves are shown in Figure [2] 
The differential comparison star measurements 
AF(Ci—C2) show a relatively large amount of 
scatter (see Fig. CO due to the faintness of the 
comparison stars relative to IRAS 17436+5003 
(V=7.0) and, from 1994—2007, the faintness of 
the initial C2 star relative to Ci. There is ev¬ 
idence of a general increase in brightness of Ci 
by 0.02 mag (V) from 1994 to 2007; this is rela¬ 
tively small and would not affect our main results. 
Inspection of the seasonal light curves of IRAS 
17436+5003 shows in general a cyclical variation 
with a period in the range of ~40—50 day. At 
times, however, it appears as though the variation 
has gone away, such as at the ends of the 2009 and 
2011 observing seasons (e.g., Fig. m- 

The ( V—Rc) color of the system shows a vari¬ 
ation of 0.07 mag peak-to-peak (except for four 
outlying points), with most of the data within a 
range of 0.05 mag; the (B— V) color range is larger, 


reaching 0.12 in 2010, with most of the data within 
a range of 0.07 mag. An inspection of the light 
and color curves shows that, in general, there is a 
trend of color with brightness, being redder when 
fainter. This can perhaps be seen more clearly in 
a plots of AV versus A (V—Rc) and AV versus 
A (B— V), as shown in Figure [3] 

4.2. IRAS 18095+2704 

The VRc light curves for IRAS 18095+2704 are 
shown in Figure |4] They display a general trend 
of increasing brightness over the two decades of 
observations, increasing by 0.30 mag in V from 
1994—2012 and 0.26 mag in Rc from 1995—2012. 
At the same time, there is very little change in 
the seasonal (V—Rc) color of the system, <—0.01 
mag, apart from what appears to be a slight zero- 
point shift of —0.01 mag that arose when we 
changed from the initial to the newer CCD in 2008. 
We assume that this difference arises from some 
unresolved imprecision at this level in the stan¬ 
dardization of the two systems for this star. 

Superimposed on this brightening trend is a 
clear cyclical variation of greatly varying ampli¬ 
tude. The amplitude (peak-to-peak) varies in V 
from 0.05 (1997) to 0.14 mag (2003) and in Rc 
from 0.02 (1997) and 0.04 (1999) to 0.11 mag 
(2003). For years in which we have good coverage 
with both filters, the average ratio of the ampli¬ 
tude in the two (Rc compared with V) is 0.86. 
In the six years in which we observed consistently 
with all three filters, the average ratio of the am¬ 
plitude in B to V is 1.29. Samples of the more 
recent B and color data, showing the cyclical vari¬ 
ations, are shown in Figure [2] 

Looking at the (V—Rc) color (Fig. [4]), it can 
be seen that in general it follows the cyclical vari¬ 
ation, being redder when fainter. This is perhaps 
seen even more clearly in the (B— V) data (Fig. [2]). 
Plotted in Figure [3] are the AF versus A (V—Rc) 
and A V versus A (B— V) data sets, with the mono¬ 
tonic trend in A V brightness first removed using 
a low-order polynomial fit. The ( V—Rc) color be¬ 
tween the two different CCD systems have been 
adjusted by the small offset of —0.01 mag men¬ 
tioned above. There is a clear linear trend, with 
the object being redder when fainter. The range in 
A (V—Rc) is 0.06 mag, with most within a range 
of 0.04 mag, the range in A (B— V) is 0.09 mag, 
with most within a range of 0.05 mag. 


4 


4.3. IRAS 19386+0155 

IRAS 19386+0155 also displays a general in¬ 
crease in brightness over the interval of observa¬ 
tion (1994—2007), as seen in Figure[5] It shows an 
increase of 0.13 mag in V over 14 years. During 
this time, there is a suggestion that the ( V—Rc) 
color was bluer by 0.02 mag in the first several 
years, but this result is uncertain since it is based 
on only a small number of Rc observations in the 
years 1995-1997 compared with the later observa¬ 
tions from 1999 and from 2002 through 2007. 

A cyclical variation is also seen in the seasonal 
data, with a varying amplitude. The amplitude 
ranges in V from 0.07 (1996, 1997, 1999) to 0.20 
mag (2003). There are many fewer Rc observa¬ 
tions in the early years, but from 2003 through 
2007, the ratio of amplitudes (Rc compared with 
V) is 0.84. 

When comparing the seasonal ( V—Rc ) color 
variations with the seasonal cyclical variations in 
brightness (Fig. 0, one can see a general pattern 
of the object being redder when fainter. This can 
be seen clearly in Figure [3l following the removal 
of the brightness trend in V using a low-order poly¬ 
nomial fit. The range in color is small, ~0.06 mag 
(ignoring three outliers). 

4.4. IRAS 19475+3119 

The light curves of IRAS 19475+3119 show no 
long-term trends in brightness, although they do 
show some small variations in the mean seasonal 
brightness and some changes in the seasonal range 
in brightness. These are seen in Figure [6] The 
seasonal range in V brightness varies from lows of 
0.07 (1995, 1997) and 0.08 (2010) mag to highs of 
0.19 (2009) and 0.17 (2003) mag. The range in Rc 
is typically less, with an average ratio (Rc com¬ 
pared with V) of 0.94, and the range in B is typi¬ 
cally more in the six years in which it is observed, 
with an average ratio (B to V) of 1.29. Visual 
inspection of the seasonal light curves shows cycli¬ 
cal variations in brightness, and in several of the 
seasons (1996, 1998, 2003, 2008, 2009, and 2010) 
one can discern a cycle length of 30—50 day (see 
Fig. HJ. 

The overall color of the system does not vary 
much, with the A (V—Rc) value mostly within 
the range of 0.035 mag and A (B— V) value mostly 
within the range of 0.04 mag. When one examines 


the seasonal light and color curves, one can see a 
trend, with the object fainter when redder (see 
Fig. 0. Such a trend can be seen clearly when 
looking at the overall brightness-color curves, as 
shown in Figure [3] 

5. COMBINED LIGHT CURVES & PE¬ 
RIOD ANALYSES 

All four of these PPNs are known to be variable 
in light and have been included in prior studies by 
other investigators. This is particularly true of 
the brightest of these sources, IRAS 17436+5003, 
which has been observed by several investigators, 
particularly Fernie. The other three sources have 
been observed extensively by Arkhipova and col¬ 
laborators, as described below. Their UBV obser¬ 
vations are carried out with a single-channel pho¬ 
toelectric photometer using a 27" aperture, and 
they tied their observations to the UBV system 
by observations of standard stars. Data are also 
available for two of the objects in the V band from 
the All Sky Automated Survey fASAS:IPo imanskil 

2002IM 

The availability of these additional observations 
affords an opportunity to enlarge the data sample, 
either in time or density or both, with the promise 
of better documented light curves. These in turn 
can yield better determined results for the period 
analysis, including the determination of multiple 
periods. This is, of course, predicated on the abil¬ 
ity to accurately combine the different data sets. 
All of the investigators made differential observa¬ 
tions, as we did, and to combine them requires ac¬ 
curate standardized measurements of the various 
comparison stars. Since different investigators use 
different comparison stars and different aperture 
sizes, this can potentially lead to slight zero-point 
offsets in the magnitudes in the different data sets. 
To correct for such offsets, we compared the vari¬ 
ous data sets where they overlapped in time, to de¬ 
termine and correct for any such offsets. These are 
discussed for each star. These individually-derived 
offsets are well established, based on a minimum 
of 18 nights, and are listed in Table 0 along with 
access to the published data sets. We briefly dis¬ 
cuss these other studies, star by star, along with 
the resulting combined light curves, and we then 
proceed to the period analyses of each. 

2 http: //www. astrouw .edu.pl/asas/ 
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Period studies were carried out for each of the 
four P PNs using the period -finding program Pe- 
riod04 ( Lenz fc Bregerll2005ll . This uses a Fourier 
transform, and it easily allows for the investigation 
of multiple periods. We began with the large, com¬ 
bined V data sets and then examined the data in 


subsets. For IRAS 18095+2704 and 19386+0155, 
the general trends were first removed. The light 
curves each show some variations in the yearly 
means. We first searched the data sets for pe¬ 
riodicity with these included, in case there were 
longer-term periodicities that these represented. 
For none of these PPNs was a significant longer- 
term periodicity found, so we then normalized, or 
adjusted, the yearly data by the yearly means and 
then did the final period analyses on the adjusted 
data sets. 


5.1. IRAS 17436+5003 

5.1.1. Combined light curves 


Due to its brightness and its unusual distinction 
of being a high-latitude object with a supergiant 
spectrum, IRAS 17436+5003 (HD 161796) became 
the target of many variability studies even before 
it was identified as a PPN. The most extensive 


of these was by Fernie an d col l aborator (|Fernje 


1983 


1986, 

1993 

, 1995ft 


19891 . 1990allbl . Il99ll : iFernie fc Seasrei 


, who carried on a long series of obser¬ 


vations including most years from 1980 through 
1998. They always include V observations and 
from 1986 to 1998 the observations are UBV. Be¬ 
ginning in 1986, the observations were made with 
an automated photoelectric telescope at a good 
site, and the number of annual observations in¬ 
creased significantly. 

We combined with our V observations the long 
and consistent data set of Fernie’s. In compar¬ 
ing the data of Fernie’s with ours during the over¬ 
lap from 1984 to 1998, we find that there are 22 
dates on which we observed on the same night, 
and from these we find an average offest of 0.000 
±0.002. Thus there is no offset to apply. We ex¬ 
te nded the coverag e sligh tly by including the data 
bv IPercv fe Welch ( 198 ill from 1979—1980; com¬ 
paring these with those of Fernie from 1980 indi¬ 
cated no offset. 


The combined Flight curve is displayed in Fig¬ 
ure [3 The agreement of the different data sets is 
good in the regions of overlap and very usefully 


increases the temporal coverage. The light curve 
shows a general increase in brightness of 0.04 mag 
from 1984 to 2012, but the increase is not strictly 
monotonic. The yearly amplitude is seen to vary 
throughout this interval by a factor of three to 
four. 


5.1.2. Period analyses 


The combination of our data with those of 
Fernie and Percy & Welch gives us a relatively- 
long time interval of 34 years of data for IRAS 
17436+5003. This affords a good opportunity to 
not only find the best period(s) for the entire data 
set, but to also examine the observations in smaller 
time intervals to investigate period changes. The 
data sets were first adjusted to the mean bright¬ 
ness of each year. Years with fewer than ten data 
points were not included. 

The analysis of the entire V data set (set A) 
revealed in several closely-spaced frequencies, as 
can be seen in Figure [ 8 ] The analysis revealed 
dominant periods of Pi = 45.15 ± 0.01 days, P 2 
= 47.36 days, and P 3 = 46.75 days. There are in 
total eight significant periods in the data; we have 
listed the first six in Table [ 7 ] The data, phased to 
Pi = 45.15 days, are also displayed in Figure [ 8 ] 
They show a reasonably good phased light curve, 
considering the observed range in amplitude seen 
in the light curve. The fit of the periods and ampli¬ 
tudes to the normalized light curve is remarkably 
good for most of the time interval, as can be seen 
in the lower five panels of Figure [ 7 ] It is only in 
the last decade, 2003—2012, that the fit in some 
seasons is not good. Some of this may be due 
to the observed complexity of the light variations 
and some may be due to changes in the periods. 
Comments about the light curve complexities and 
period changes were frequently made by Fernie. 
Fernie found dominant periods ranging from 38 to 
62 days when he examined the observations in one 
or two year intervals, although most were in the 
range of 40—46 days. Fernie also drew attention 
to so me epochs in which the star appeared to not 
vary ( Fernie![Ti)9.'ji ). Some of this complexities seem 
to be resolved by our inclusion of multiple periods. 

The data were then examined in subsets and 
the results are also listed in Table [3 The anal¬ 
ysis of the Fernie and Percy & Welch data (set 
B, 1979—1998) results in Pi = 45.08 ± 0.01 days, 
P 2 = 46.54 days, and P 3 = 47.45 days; these are 
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nearly the same as the first three periods seen in 
the entire data set. This is perhaps not surprising, 
since their data comprise two-thirds of the number 
in the entire data set. The analysis of the com¬ 
plete VUO data (set C, 1994—2012) results in Pi 
= 47.30 ± 0.02 days, which agrees with P2 in the 
entire data set (A), and a value for P 2 that agrees 
with P4 in the entire data set. Similar results were 
found in the VUO Rc data (set D, 1997—2012). 
We also looked at subsets of the VUO V data in 
five or six year intervals. The 1994—1999 data (set 
E) has Pi = 44.09 and P 2 = 42.38 days; the sec¬ 
ond of these agrees with P3 of the entire VUO V 
data set (C). The 2003—2007 data (set F) has Pi 
= 41.47 and P 2 = 56.87 days; this value of Pi 
does not appears in any of the other data sets. 
The 2008—2012 data (set G) has Pi = 47.42 and 
P 2 = 57.70, and P3 = 38.14 days; Pi agrees with 
the value found for Pi in the complete VUO data 
set (C) and P 2 is similar to the value of P 2 found 
in the 2003—2007 data set. 

Thus, while there are several common periods 
found in many of the data sets, particularly 45.1 
and 47.4 days, there are also values of P that are 
found in only certain time intervals (sets E, F). 
This suggests that there are indeed changes in the 
dominant periods with time. However, there does 
not appear to be a monotonic increase or decrease 
in the period over time, as seen by the compari¬ 
son of the Pi values in subsets E, F, and G. One 
can also meaningfully compare the ratio of the 
two dominant periods, particularly in the longer 
data sets (A, B, C, D). They are all close in value, 
ranging from 1.03 to 1.06, with an average value 
of 1.05. 


5.2. IRAS 18095+2704 


5.2.1. Combined light curves 

IRAS 18095+27 04 was observed bv lArkhioova et a! 
di 993I. I2000L l20inh 'or 19 years, beginning with 
only a few observations in 1990—1992 but then 
more frequently through 2008. They quote an un¬ 
certainty of ±0.01 mag in B and V. Their observa¬ 
tions show a clear trend of increasing brightness, 
as we found, with an increase of 0.37 mag in V. 
0.35 mag in B, and 0.32 mag in U over 19 years, 
with the (B—V) color staying approximately con¬ 
stant. Examining some older photographic ob¬ 


servations, they find that this brightness trend 
extends back to the earliest measurements from 
the mid-1930s. From their 2000—2008 observa¬ 
tions, they find a period of 109+2 days with a 
semi-amplitude of 0.03 mag, and they note that 
during this cyclical variation the system is bluer 
when brighter. 

To combine our V data with theirs, we first 
compared the two data sets. There are 24 nights 
in common between the two data sets, with a sys¬ 
tematic difference (VUO — Arkhipova) of —0.054 
mag and another 50 nights in which observations 
from the two data sets were made on successive 
nights, with a difference of —0.042. This yielded a 
weighted value of —0.048 mag, with those on suc¬ 
cessive nights receiving half weight. Possible rea¬ 
sons for this offset include uncertainties in stan¬ 
dardization of the measurements, including the 
comparison star, and perhaps, due to their larger 
aperture, an unseen star in their comparison star 
aperture. 

As part of the ASAS sky survey, 69 observations 
of IRAS 18095+2704 with a typical uncertainty of 
±0.03 mag were made during 2003 and part of 
2004. The 2003 data show a well-delineated cycli¬ 
cal pattern over almost two cycles with a period 
of ~105 days. However, we did not combine these 
data with ours. 

A combined V light curve was formed by adding 
the Arkhipova et al. data set from 1993 to 2008 
to ours, using the offset determined. This is dis¬ 
played in the top panel of Figure [9l and shows an 
increase in brightness of 0.40 mag over this 20- 
year interval. To prepare the combined data set 
for analysis, we first removed the trend in each 
data set separately using a low-order polynomial, 
and then adjusted the seasonal light curves to the 
same mean levels. This adjusted V light curve is 
shown in the bottom panel of Figure |9] Close ex¬ 
amination of the seasonal light curves shows good 
agreement between the two data sets and serves 
to better document the cyclical variability. 

5.2.2. Period analyses 

We began with a period analysis of the com¬ 
bined, adjusted Flight curve of IRAS 18095+2704, 
which covers the longest time interval, 1993 to 
2012, and has the largest number of data points 
(set A). The frequency spectrum showed several 


7 













relatively strong peaks, and the subsequent anal¬ 
ysis indicated multiple periods in the data. The 
strongest period was Pi = 113.3 ± 0.1 days, with 
additional periods of 97.9, 103.0, and 108.4 days. 
The frequency power spectrum and the phase plot 
with Pi are shown Figure 0 and the presence of 
these additional periods is seen in the frequency 
spectrum. Six significant periods are found in the 
data, and these are listed in Table 0 

In Figure [U] is also shown the fit of the first 
four periods, along with their associated ampli¬ 
tudes and phases, to the observed light curve. The 
fit is reasonably good, given the complex nature 
of the light curve and the long interval of obser¬ 
vations (20 years), and improves slightly with the 
inclusion of the two additional significant periods. 

Subsets of this V light curve and of some of 
the light curves with other filters were also investi¬ 
gated. We began by dividing the combined V data 
set into two ten-year time intervals. For the time 
interval 1993—2002 (set B), the analysis reveals 
Pi = 98.2 days, close to P 2 in the full V data set, 
and an additional significant period of 107.5 days 
(close to P4 in the full data set). For the interval 
2003—2012 (set C), there is a strong period at Pi 
= 111.8 days, close to Pi in the full V data set, and 
additional periods of 110.2, 103.1 (close to P 3 in 
the full data set), 68.3, and 97.2 days; the periods 
of 111.8 and 110.2 in the 2003—2012 V data set 
are close to and beat against one another. Note 
that a formal analysis of the ASAS 2003—2004 V 
light curve results in a period of 110 ±1 days. It 
can be seen that the same periods are not dom¬ 
inant throughout the 20-year observing interval, 
although almost all of those in the two 10-year in¬ 
tervals are found in the combined data set. Next 
we investigated the B light curve, formed by re¬ 
moving the trends from the Arkhipova et al. and 
the VUO data, including the determined B offset, 
and adjusting to the individual seasons. Excluded 
were two years with fewer than 10 observations. 
This data set (D) yielded periods in agreement 
with those found for the combined V light curve, 
Pi = 113.5 days and additional periods of 103.0, 
154.5, 97.9, and 116.3 days, and with larger ampli¬ 
tudes than found for the V analysis. We also inves¬ 
tigated the VUO Rc light curves from 1998—2012 
(set E) and 2003—2012 (set F). Some of the same 
periods appear, but not exactly, and the Rc light 
curves have smaller amplitudes. These periods are 


also recorded in Table [7] Based on the combined, 
full Flight curve (set A), we find a ratio of P2/Pi 
= 0 . 86 . 

5.3. IRAS 19386+0155 

5.3.1. Combined light curves 

2010h also ob¬ 
served IRAS 19386+0155 for 19 years, again be¬ 
ginning with only a few observations in 1990—1992 
but then more frequently through 2008. They 
quote an uncertainty of ±0.01 mag in B andV. 
Their observations show a clear trend of increas¬ 
ing brightness with time. They determined an in¬ 
crease in V of 0.2 mag but find that surprisingly 
the (B—V) color gets redder by 0.1 mag during 
this time. From the 2000—2008 observations, they 
derived a dominant period of 102 day and addi¬ 
tional periods of 98 and 106, yielding a period ratio 
P2/P1 = 0.96. They note that during this cyclical 
variation the system is bluer when brighter, as we 
found. 

To combine our V data, we compared the two 
data sets and found that there are 12 nights in 
common, with a systematic difference (VUO — 
Arkhipova) of +0.060 mag and another 26 suc¬ 
cessive nights of observations, with a difference of 
+0.061. This yielded a weighted value of +0.060 
mag, with those on successive nights receiving half 
weight. 

IRAS 19386+0155 was also observed with in the 
ASAS sky survey from 2002 through 2009. We ex¬ 
amined the data set of 766 V observations, which 
had a average uncertainty of ±0.035 mag. Com¬ 
paring their data with ours on similar nights lead 
to an offset of —0.055 mag. Given the significantly 
larger uncertainty in the data, we chose not to 
combine them with our data; we did, however, in¬ 
vestigate them separately for periodicity. 

A combined V light curve was thus formed by 
adding the Arkhipova et al. data from 1993—2008 
with ours, with the offset included. This is shown 
in Figure [lUl The data sets show good agreement, 
both in the overall trend in brightness and in the 
seasonal cyclical variations present. 

5.3.2. Period analyses 

The combined VUO and Arkhipova et al. V 
light curve of IRAS 19386+0155 was analyzed for 
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variability. The trend of increasing brightness was 
first removed by fitting it with a low-order poly¬ 
nomial. The light curve analysis shows two dom¬ 
inant periods in the data, with Pi = 101.8 ±0.1 
days and P2 = 98.6 ±0.1 days. The frequency 
spectrum and phased light curve for Pi are shown 
in Figure E] These are th e same periods found 
by [Arkhipova et al.1 ; 2010 i based on their data 
set alone. Three additional significant periods are 
found in the combined data, and these are listed 
in Table 0 Similar results were found, whether 
or not we adjust the data to the means of the in¬ 
dividual seasons. We have listed in the table and 
shown in the figures the results with data adjusted 
(set A). The fit is reasonably good over this 16- 
year interval. Based on the combined, full V light 
curve (set A), we find a ratio of P2/P1 = 0.97. 
We then investigated the periodicities using sub¬ 
sets of the adjusted V data — from 1993—2000 (set 
B) and 2001—2008 (set C) — and also the VUO 
Rc data from 2002—2007. They result in domi¬ 
nant periods ranging from 99 to 103 days, similar 
to those found in the combined, full V data set. 
We also investigated the ASAS V data set from 
2002—2009, which yielded a dominant period of 
102.3 ± 0.1 days, with secondary periods of 83 and 
101 days. The 2002—2005 ASAS data, which ap¬ 
pear the most cyclical, yield a dominant period of 
102.8 ± 0.2 days. These ASAS r esults agree with 
both the period value cited by Arkhipova et al 


( 2010l) of 102±1 day based upon the 2003—2004 
ASAS data alone, and they agree with the values 
found from our larger and more precise data sets. 


5.4. IRAS 19475±3119 

5-4-1- Combined light curves 

IRA S 19475±3119 was observed bv l Arkhipova et~al 
( 2006 ) over a four-year interval from 2002 to 2005, 
resulting in 104 UBV measurements. They re¬ 
ported that the object showed semi-regular bright¬ 
ness variations with a maximum range of 0.17 mag 
in V. and a larger range in the shorter bandpasses. 
This agrees with the variations that we find. They 
also reported that they did not find a statistically- 
significant period for the entire data set, but pe¬ 
riods were found for individual years of 43±1 day 
for 2003 and 45±1 day for 2004. Upon download¬ 
ing their archived data, we found that it included 
some more observations from the following season, 


for a total of 132 UBV measurements. 

Since their data set falls within the time range 
of ours, including their data will not extend the 
base line but will increase the density of obser¬ 
vations during that five-year interval of overlap. 
Comparing V observations from the two data sets 
revealed an offset in observations made on the 
same day of (VUO—Arkhipova) of —0.094 mag (6 
nights) and an offset of observations made on suc¬ 
cessive days of —0.087 mag (12 nights). This yields 
a weighted offset of —0.091 mag, where those made 
on successive dates received only half weight. The 
combined V light curve with the offset included is 
displayed in Figure HU Visual inspection season 
by season shows good agreement between the two 
data sets. 


5.4-2. Period analyses 


A period analysis of IR AS 19475±3119 was first 
carried out on the combined Flight curve. We an¬ 
alyzed the light curve both with and without ad¬ 
justing the data to the seasonal means. The same 
results were obtained in both cases, with two sig¬ 
nificant periods, Pi = 40.95 ±0.01 days and P2 = 
38.85 ±0.01 days (set A). An analysis of the nor¬ 
malized VUO Rc light curve (set B) also resulted 
in two similar periods. These period results are 
listed in Table 0 The frequency spectrum and 
the phased light curve plot based on Pi for the 
combined and adjusted V data are shown in Fig¬ 
ure [8] Both the V and R light curves also gave 
evidence of a third period, P 3 = 48.1 days, with a 
S/N = 3.9, slightly below our significance level of 
4.0. However, given that it is seen in both of the 
data sets, we regard it as significant. These each 
yield a ratio of P 2 /P 1 = 0.95. 

In Figure fill is displayed the adjusted Flight 
curve for the combined data, fitted with the three 
periods and associated amplitudes and phases. 
The fit for this star is not as good as that found 
for the other three, and we attribute much of this 
to the complexity of the pulsations of this star. 
A better fit could likely be obtained with time- 
dependant amplitudes. PPNs h ave been shown to 
produce shocks as they pulsate ( Lebre et al. I ll996t 
Zacs et adll2009h . A spectroscopic and photomet¬ 
ric study of the similar (F5 I, P~39 days) but 
carbon-rich PPN IRAS 07134±1005 (HD 56126, 
CY CMi), with a similarly complex light curve, 
displayed strong evidence for shock waves in its 
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atmosphere ( Barthes et al . ISqo3). 

We also analyzed the light curves in subsets and 
using the other filters. Analyses of the V (set C) 
and Re (set D) light curves from 2002—2012 yield 
similar periods Pi and P2 to the full data sets. 
Only the Plight curve for the interval 1994—1999 
(set E), based on many fewer points, gave signifi¬ 
cantly different periods of 36.7 and 36.1 days, val¬ 
ues not found in the other data sets. This suggests 
a change in the period from the earlier to later 
time intervals. However, this result is based on a 
smaller data set over fewer years, and thus these 
periods are less secure than the periods found in 
the other data sets for this object. 

6. DISCUSSION 

We have analyzed light curves that span time 
intervals from 16 to 34 years for four bright, O-rich 
PPNs, all with F spectral types, Dominant periods 
were determined for all four, with values ranging 
from 41 to 113 days. Similar, although less robust 
values, had been found in previous studies. These 
values fall within the range found for other F—G 
spectral type PP Ns, in particu l ar the 12 C-rich 
ones analyzed bv iHrivnak et all ( 2010h that span 
the period range of 38 to 153 days. A similar range 
was fou nd for eight C-rich PP Ns in the Magellanic 
Clouds ( Hrivnak et alJl2015al Paper III). 

We also analyzed the light curves within subsets 
of five to ten year time intervals. These analyses 
often resulted in different dominant period values 
in different time intervals, indicating that signif¬ 
icant period changes do occur. In a majority of 
cases, Pi increased with time, but not in all cases. 
In a previous period analysis of two C-rich PPNs 
with similarly long sets of observations, we found 
that in one case there was no signific ant cha nge in 
Pi and in the other Pi decreased (IHrivnak et al 
20131 1. With this present small sample, we con¬ 


clude that no clear trends in period change emerge 
when comparing the data sets for an individual 
star in different time intervals. The results of 
the larger sample of 12 PPNs had suggested that 
Pi decreases with time ( Hrivnak et alJl201flh . but 
that result is not seen over the longer two-decade 
time span of these studies. 

Multiple periods were found for all four PPNs, 
and they were always close in value to the dom¬ 
inant period. Similar results have been found 


i n detailed light curve stu dies of other PPNs 
( Hrivnak et al. 2013 . 20 1 5 all . The ratios of P2/P1 
for these four range from 0.86 to 1.05, with an 
average value of 0.96. These are summarized in 
Table [HJ They are s imilar to the r a tios fo und in 
previous studies by IHrivnak et al. (120131 Paper 
II) for two PPNs ('0.95b bv lHrivnak et all ( 2015al 
Paper III) for s e ven PPNs ( 0.86—1.16), and by 
lArkhipova et al.l ( 2010 . 12011 ) for several PPNs 
and other post-AGB stars (0.93—1.06). Also, sim- 
ila r values we r e foun d for other post-AGB objects 
bv llviss et al.l (12007 1 for three objects (0.90—0.94) 
and by Van Winckel et al. (2009) for two others 
(0.91, 0.96). Thus these seem to be typical pe¬ 
riod ratios for these low- and intermediate-mass 
post-AGB objects. They presumably relate to the 
pulsational modes excited, although these have 
not been identified. 

Pulsation periods and amplitudes, together 
with the determined temperatures, can potentially 
be used with stellar pulsational models to deter¬ 
mine th e masses and luminosities of these evolved 
objects ( Fokin et al. 2001 : Aikawal 2010l h These 
derived values can be compared to theoretical 
mass—lumino sity relationships de rived f or post- 


AGB stars (Blocker Il995l : Vassiliadis & Wood 


19941 1. At present, there are only a few investi¬ 


gations of pulsation of post-AGB stars, and t he 
agreem ent with observations is not good. lAikawal 
( 2010 1 ran linear, radial pulsation models of post- 
AGB stars with masses of 0.6 and 0.8 M 0 and 
presented graphically the results for the temper¬ 
ature range 5000—7000 K and log g in the range 
of 0.0—2.0. We examined the results for T e s = 
6300 K and T e s = 7000 K and log g in the range 
of 0.2 to 1.0, the closest parameters to our stars. 
The fundamental modes were all stable ex cept f or 
those with T e g = 6300 K and log g < 0.4. lAikawa 
( 2010 ) also carried out on-linear analyses of the 
unstable modes. However, the calculated values 
had periods that were too short (days to a few 
weeks) and amplitudes that were too small by a 
factor of three when co mpared to our ob served 


results. The models by Fokin et al. ; 2 001, : 


were 

non-linear models of masses of 0.6 and 0.8 M 0 
but of temperatures in the range of 5600 to 6000 
K, cooler than our stars. We previously compared 
them with our observational results of two cooler 
PPNs, but they also resulted in periods that were 
too short and pulsational amplitudes that were 
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too small ( Hrivnak et al . H20H). This commends 
further pulsational modeling of post-AGB stars in 
the temperature range of 6000—8000 K, . 

One can also make other comparisons with the 
results of our study of 12 C-rich PPNs in the Milky 
Wav Galaxy. For th at sample of F—G stars, we 
( Hrivnak et al . lEmJ found a monotonic, approxi¬ 
mately linear decrease in period (P) with effective 
temperature (T e //) over the range of 5000 to 8000 
K. While we do not have a long enough observing 
interval to see an individual star change in T e ff 
or P, we can interpret this trend of decreasing 
P with increasing P e // to represent the average 
expected evolution of these stars as they evolve 
from large, cool AGB stars to the small, hot cen¬ 
tral stars of PNs. Since they evolve toward higher 
temperatures at approximately constant luminos¬ 
ity, they must be decreasing in size and increasing 
in density (ongoing mass loss is small). There¬ 
fore, assuming radial pulsations, we would expect 
a decrease in period as they evolve. A compar¬ 
ison with five C-rich PPNs in the Large Magel¬ 
lanic Cloud (LMC) showed that they also follow a 
similar monotonic trend, although perhaps offset 
slightly to lower T e ff or shorter P (observations 
of more periodic LMC PPNs are needed to con¬ 
firm this offset). In Figure IT2l is plotted the P and 
T e ff values for the four O-rich PPNs examined in 
this study, together with the 12 C-rich ones from 
the Milky Way Galaxy. Note that we have in¬ 
cluded two points each for IRAS 17436+5003 and 
IRAS 19475+3119, since there are two different 
T e ff values published for each. A straight line fit 
to the data points for these 4 O-rich PPNs would 
be nearly vertical and possess a much steeper than 
that found for the C-rich ones. However, this 
present sample is small and all of the points, ex¬ 
cept one of those for IRAS 17436+5003, fall within 
the range of the C-rich data points. When we com¬ 
pare the pulsation amplitude (AV) with their pe¬ 
riods and temperatures (see Table 0, they agree 
well with the monotonic trends found for the C- 
rich sample; all of the objects with periods shorter 
than 120 days and temperatures higher than 6000 
K have maximum amplitudes <0.22 mag. This is 
shown in Figure flTTl 

As noted above, the P—T e ff relationship for 
these four O-rich PPNs appears to have a much 
steeper slope than is found for the 12 C-rich PPNs. 
This may not be significant, since it is based 


on only four objects over a smaller temperature 
range (1000 K), and the O-rich PPNs generally fall 
within the distribution of the C-rich sample. If it is 
real, one suggestion is that it relates to the differ¬ 
ence in metallicity between the two longer-period 
objects (P«110 day, [Fe/H]«—1.0) and the two 
shorter-period objects (P?»45 day, [Fe/H]«—0.25). 
This would indicate that the longer-period ones 
derive from a more metal-poor population and 
thus are older and of lower m ass. Th e non -l inear 
radiative pulsation models bv iFokin et al.1 ( 2001 ) 
indicate longer periods for lower-mass stars, which 
is consistent with what we find here. However, this 
is only suggestive and may not be relevant since 
the models were for somewhat cooler (5600—6000 
K) stars than these four. If this difference in slope 
is confirmed by a larger study, then further theo¬ 
retical models to investigate this would be desir¬ 
able. 

Two of the four objects in this sample show 
significant long-term trends in brightness over the 
observing interval, both getting brighter: IRAS 
18095+2704 by 0.30 mag over 19 years or 0.016 
mag yr -1 and IRAS 19386+0155 by 0.13 mag 
over 14 years or 0.009 mag yr” 1 . Similar mono¬ 
tonic trends, some decreasing and some increas¬ 
ing, have been found in the light curves of other 
PPNs of F—G spectral types. Including these four, 
periods have been determined for a total of 26 
PPNs of F—G spectral types, 19 in the MWG 
and 7 in the Magellanic Clouds, for which ob- 
servations cover intervals of at least nine years 
j Hriynak_et_a 1 2010: lArkhinova et al. 2010L 2011 


Hrivnak et al.l2015aj ). Of these, five show mono¬ 

tonic increases of ^0.01—0.02 mag yr -1 , two show 
monotonic decreases of ~0.01 mag yr -1 , and one 
shows a light curve that is constant in brightness 
for six years, then shows a sudden drop of 0.12 
mag in one year, and then a gradual increase of 
0.019 mag over nine years. Two others show some 
larger changes (~0.1 mag) that are not monotonic. 
The majority (16 of 26) do not show much change 
in the seasonal brightness levels. 

Gradual changes in the brightness could be due 
to secular changes in the dust opacity. PPNs are 
surrounded by an expanding envelope of circum- 
stellar gas and dust, which formed from mass 
loss on the AGB. This circumstellar envelope 
can both dim and redden the emerging light. 
The sudden drop and reddening witnessed in one 
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PPN (IRAS 19500—1709) has been attributed 
to the for mation or transit o f a circumstellar 
dus t cloud ( Hrivnak et, al ■ 2010 : lArkhioova et al 


2010). Both IRAS 18095+2704 and 19386+0155 


are significantly redde ned (see Table [TJ (B—V) 
= 0.23-0.32 for F2-5I. Icralkoooh . However, the 
lack of color change in these two PPNs when they 
increase in brightness would require that the dust 
opacity along the line of site be due to large (> 1 
/im) grains, which wo uld appear grey and not color 
selective, as noted bv I Arkhipova et all ( 2010h . 


All four of these PPNs are redder, and thus 
cooler, when fainter. This is what is com- 
monly found in PPNs of F—G spectral types 


(Arkhipova et al. 2010i 

2011: 

Hrivnak et al. 2010L 

2015a). Hrivnak et al. 

(2013) found in the de- 


tailed study of two PPNs that included contem¬ 
poraneous light, color, and velocity curves, that 
the stars were smallest when brightest and hottest. 
This differs from what is found in Cepheid vari¬ 
ables, in which there is a phase lag of ^0.25 be¬ 
tween the maximum in light and the minimum of 
size. This temperature and size correlation will 
be investigated for three of these O-rich PPNs 
through a coordinated light, color, and velocity 
study, for which observations are presently being 
made. 


The specific changes in color can be correlated 
to specific changes in temperature for these four 
objects. The maximum ranges in colors are listed 
in Table [8] Using the observed (B— V) changes 
and the temperatures determined from the spec¬ 
troscopic analyses, these color changes were trans¬ 
formed to temperat ure chan ges using the color- 
temperature table of Coxl ( 21)1)1) . Table 15.7). This 
resulted in maximum temperature changes of 400 
to 770 K, which are attributed to pulsation. Of 
the two PPNs with multiple spectroscopically- 
determined temperature values, for one of them 
(IRAS 17436+5003) the difference between the 
two spectroscopic temperatures is less than the 
range determined from the color range, and for 
one of them (IRAS 19475+3119), the difference 
is slightly greater, but within the uncertainty of 
the temperature determinations. We investigated 
the epochs of the high-resolution spectra used to 
determine the temperatures, to see if they could 
be correlated with specific phases in the light 
curves. For IRAS 17436+5003, the higher tem¬ 
perature spectrum was obtained near maximum 


brightness and bluest color (see Fig. [7]), while the 
lower temperature spectrum was obtained when 
there were no contemporaneous light curves. For 
IRAS 18095+2704, the spectrum was taken near 
maximum brightness and bluest coloJl (see Fig. 0), 
and for IRAS 19386+0155, between maximum and 
average brightness (see Fig. Col). Of the two spec¬ 
tra of IRAS 19475+3119, the higher tempera¬ 
ture one was obtained near maximum brightness 
and the lower temperature one at an epoch with 
few photometric observations and no clear phas¬ 
ing of the brightness. Thus, for the two objects 
with multiple spectroscopic temperature determi¬ 
nations, the higher temperatures are associated 
with the greater brightness in the light curves, con¬ 
sistent with the brightness-color correlations found 
from the light curves. 

7. SUMMARY AND CONCLUSIONS 

In this paper, we have carried out a detailed 
light curve and period study of four bright PPNs. 
They each had previous light curve and period 
studies, but by combining these data sets with 
ours, we have significantly increased the sample 
size and the time intervals of observation. The pri¬ 
mary results are listed below. They serve to help 
elucidate the pulsation properties of intermediate- 
and low-mass post-AGB during this short transi¬ 
tion stage from the AGB to the PN phases. 

1. The four vary in light with changing am¬ 
plitudes; these amplitudes vary between 0.05 and 
0.25 mag, peak to peak. 

2. Two of them also show long-term trends 
of increasing brightness, with values of 0.13 and 
0.30 mag over 14 and 19 years, respectively, cor¬ 
responding to 0.01—0.02 mag yr . Long-term 
trends have been observed in about 30% of the 
well-studied PPNs, with five showing clear mono¬ 
tonic trends of increased brightness and two of de¬ 
creased brightness. We attribute these to changes 
in the circumstellar dust opacity. 

3. A dominant period was found in each of 
the four, with values ranging from 41 to 113 days. 
Similar periods had previously been found for 
three of these based in smaller data sets. These 

ISahin et al.l 120111 ') obtained two separate spectra, observed 
10 months apart. Since they based the analysis primarily 
on the first of these, this is the epoch we investigated. 
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values fall in the period range found previously 
for PPNs of F—G spectral types of 38 to 153 days. 

4. In each of the four PPNs, multiple peri¬ 
ods were determined. These account for much 
of the variation in amplitude seen in the light 
curves. In all cases, the secondary period is close 
to the primary period, with an average value of 
P 2 /P 1 of 0.96. This is similar to the range of val¬ 
ues of 0.9—1.1 found in previously studied PPNs 
with suitably long temporal data sets. These mul¬ 
tiple periods help to explain the changing light 
curve amplitudes as they beat together. They also 
help to explain the seasonally changing periods at¬ 
tributed to IRAS 17436, which instead are due to 
multiple periods affecting the light curve shape. 

5. For these four, one finds a general trend 
of the cooler ones having longer periods. How¬ 
ever, this is not a strong trend and the sample is 
small. All four of these are O-rich. A strong trend 
has been seen in C-rich F—G spectral type PPNs. 
These four O-rich ones have a steeper slope to their 
trend, but the sample is too small to determine 
if the slopes differ between PPNs with different 
chemistries. 

6. All four are redder when fainter, in agree¬ 
ment with the brightness-color relationship found 
in other well-studied PPNs. The maximum 
(B—V) color changes range from 0.055 to 0.105 
mag, which imply changes in effective tempera¬ 
ture of 400-770 K. 

We have begun photometric studies of addi¬ 
tional, fainter 0-ricli PPNs, which will help to 
enlarge the sample and increase the range in tem¬ 
perature. These should help to better define the 
period-temperature-amplitude relationships for 
0-riclr PPNs and allow a more robust compari¬ 
son with C-rich PPNs. 
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Table 1 

Program Objects 


IRAS ID 

V a 

(mag) 

(V—R c ) a ’ 6 
(mag) 

Sp.T. 

T e ff 

(K) 

log g 

[Fe/H] 

[C/O] 

Ref. c 

Other ID 

17436+5003 

7.0 

0.2 

F3 Ib d 

6600 

0.0 

-0.3 

-0.2 

1 

HD 161796, V814 Her 


7.0 

0.2 

F3 Ib d 

7100 

0.5 

-0.25 

-0.3 

2 


18095+2704 

10.0 

0.7 

F3 lb 

6500 

0.5 

-0.9 

-0.37 

3 

V887 Her 

19386+0155 

11.1 

0.8 

F5 lb 

6800 

1.4 

-1.1 

-0.35 

4 

VI648 Aql 

19475+3119 

9.3 

0.4 

F3 lb 

7750 

1.0 

-0.24 

-0.39 

5 

HD 331319, V2513 Cyg 


9.3 

0.4 

F3 lb 

7200 

0.5 

-0.25 

-0.5 

2 



a Variable. 

b Includes circumstellar and interstellar reddening. 

c Re feren ces for the spectrosc opic analyses: (1) [L uck et al.l dl990t) , (2) [Klochkov a et all d2002th (3) ISahin et al.l 
d201llh (4) |Pereira et alj d2004T). (5 ) I Are llano Ferro et al.1 d200ll~ Note that there is an earlier spectroscopic study 
of IRAS 18095+2704 bv lKlochkov^ (|1995l h but the one cited is at ~2.5 times higher resolution. 

d Classified as A7 I by I Suarez et al.l d20061) . 
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Table 2 

Standard Magnitudes of Program and Comparison Stars 1 


Object 

GSC ID 

V 

(B-V) 

(V-R C ) 

(Rc-lc) 

Run b 

IRAS 

17436+5003 

03518-00402 

7.00 


0.25 


1 




6.94 

0.35 

0.22 


2 

Ci° 


03518-00926 

10.32 

0.81 

0.46 



c 2 


03518-01037 d 

9.65 

1.05 

0.58 



c 3 


03518-01174° 

9.64 

1.03 

0.59 



IRAS 

18095+2704 

02100-00044 

10.30 


0.72 


1 




10.00 

1.02 

0.69 


2 

Ci 


02100-00387 f 

12.28 

0.55 

0.33 



C 2 g 


02100-00238 h 

12.41 

0.88 

0.54 



c 3 


02100-00003 

13.57 

0.98 

0.56 



IRAS 

19386+0155 

00483-00956 

11.09 

1.16 

0.77 

0.86 

3 

Ci 


00483-01010 

13.04 

0.81 

0.46 

0.49 


C 2 


00483-01014 

13.44 

0.94 1 

0.55 

0.56 


c 3 


00483-01082 

13.54 

0.83 1 

0.46 1 

0.51 1 


IRAS 

19475+3119 

02669-01757 

9.27 


0.41 


1 




9.32 

0.63 

0.40 


2 

C+ 


02669-02709 

11.02 

0.14 

0.08 



c 2 


02669-04063 

11.92 

1.76 

1.01 



c 3 k 


02669-04215 

12.39 

1.32 

0.75 




Uncertainties in the observations are as follows — V: ±0.01, (B—V): ±0.025, (V— Rc) : 
±0.015 mag, and (Rq—Ic)'- ±0.010 mag, except as noted. 

Standardized observations were made at the VUO on UT (1) 22 August 1995 and (2) 
26 June 2012, and are averaged together for the comparison stars, and at (3) the Kitt Peak 
National Observatory on 23 June 1994. 

c Evidence of a general increase in brightness from 1994—2007 of 0.02 mag (V). 
d HD 234482 
e HD 234480 
f TYC 2100-387-1 

Suggestion of a gradual monotonic increase in brightness of 0.025 mag ( V) from 1994 to 
2007 and a gradual increase and then decrease over a range of 0.03 mag ( V) from 2008 to 
2012 . 

h TYC 2100-238-1 
Uncertainty of ±0.02 mag. 

Suggestion of variation < 0.03 mag peak to peak. 

k Suggestion of long-term (> 8 year) cyclical variation ~0.03 mag peak to peak. 
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Table 3 

Statistics of the VUO PPN Light and Color Curves 


IRAS ID 

Years 

Number of Observations 


Average Uncertainty (mag) a 

AV 

AR C 

A(V-Rc) 

AB 

cr(AV) 

cr(ARc) 

<t(A(V-Ro)) 

<t(AB) 

17436+5003 

1994-2012 

563 

542 

450 

284 

0.006 

0.006 

0.009 

0.008 

18095+2704 

1994-2012 

530 

473 

399 

246 

0.005 

0.005 

0.007 

0.005 

19386+0155 

1994-2007 

177 

135 

109 


0.009 

0.006 

0.010 


19475+3119 

1994-2012 

481 

425 

360 

236 

0.003 

0.004 

0.005 

0.002 


The average statistical uncertainty in a single differential measurement. 



Table 4 

Differential Standard VRc Magnitudes and Colors for the Four PPNs from 1994-2007 1 


Program Star 

HJD - 2,400,000 

AV 

HJD - 2,400,000 

ar c 

HJD - 2,400,000 

A(V-Rc) 

IRAS 17436+5003 

49557.6197 

-3.334 

49873.6310 

-3.121 

49873.6309 

-0.222 

IRAS 17436+5003 

49565.5887 

-3.311 

49882.6152 

-3.085 

49882.6117 

-0.200 

IRAS 17436+5003 

49571.5873 

-3.268 

49943.6121 

-3.092 

49943.6121 

-0.245 

IRAS 17436+5003 

49615.5535 

-3.302 

49951.5862 

-3.115 

49951.5869 

-0.221 

IRAS 17436+5003 

49624.5247 

-3.256 

50263.6893 

-3.146 

50263.6889 

-0.227 

IRAS 17436+5003 

49632.5300 

-3.279 

50276.6724 

-3.103 

50276.6724 

-0.214 

IRAS 17436+5003 

49638.5021 

-3.295 

50288.6915 

-3.092 

50288.6914 

-0.201 

IRAS 17436+5003 

49643.5242 

-3.341 

50385.4892 

-3.122 

50385.4892 

-0.209 

IRAS 17436+5003 

49646.5125 

-3.363 

50591.6397 

-3.058 

50591.6431 

-0.207 

IRAS 17436+5003 

49653.4942 

-3.339 

50596.7252 

-3.074 

50596.7250 

-0.216 


Note. —Average uncertainties in the brightness are as follows: ±0.006 (AV), ±0.005 (ARc), ±0.008 mag 
(A(V—Rc)) 5 with maximum uncertainties ~2 times as great. 

a Table [4] published in its entirety in the electronic edition of the Astronomical Journal. A portion of Table [4] is 
shown here for guidance regarding form and content. 

b The 2007 observations that contain all three BVRc filters are listed instead in Table [5] 






Table 5 

Differential Standard BVRc Magnitudes for the Four PPNs from 2007-2012 a ’ & 


Program Star 

HJD - 2,400,000° 

AB 

AV 

AR C 

HJD - 2,400,000° 

AB 

AV 

AR C 

IRAS 17436+5003 

54283.6210 

-3.761 

-3.304 

-3.098 

55415.6360 


-3.301 


IRAS 17436+5003 

54287.6152 

-3.771 

-3.311 

-3.086 

55417.6187 

-3.719 

-3.287 

-3.083 

IRAS 17436+5003 

54288.6223 

-3.737 

-3.308 

-3.104 

55421.5619 

-3.714 

-3.279 

-3.077 

IRAS 17436+5003 

54301.6111 

-3.838 

-3.344 

-3.135 

55428.5967 

-3.721 

-3.298 

-3.080 

IRAS 17436+5003 

54302.6374 

-3.834 

-3.364 

-3.161 

55436.6151 

-3.802 

-3.339 

-3.116 

IRAS 17436+5003 

54304.6510 

-3.851 


-3.128 

55438.6311 

-3.803 

-3.367 

-3.137 

IRAS 17436+5003 

54304.6523 

-3.851 


-3.128 

55440.6121 

-3.836 

-3.390 


IRAS 17436+5003 

54305.6554 

-3.853 

-3.368 

-3.148 

55447.5513 

-3.875 

-3.404 

-3.153 

IRAS 17436+5003 

54311.6534 

-3.796 

-3.339 

-3.122 

55449.6145 

-3.870 

-3.394 

-3.158 

IRAS 17436+5003 

54312.5908 

-3.801 

-3.361 

-3.111 

55452.5671 

-3.862 

-3.388 

-3.151 


Note. —Average uncertainties in the brightness are as follows: ±0.005 (AB), ±0.005 (AV), and ±0.004 (ARc), with maxi¬ 
mum uncertainties ~2 times as great. The uncertainties for IRAS 17436+5003 are higher than the average and those for IRAS 
19475+3119 are lower. 

a Tablel5lis published in its entirety in the electronic edition of the Astronomical Journal. A portion of Table f5l is shown here 
for guidance regarding form and content. 

b The 2007 observations that contain only VRc filters are listed instead in Table f~4| 

c The time represents the mid-time of the V observations. The times for the B and Rq observations differ from this by 
approximately the following: IRAS 17436+5003: +0.0011 and —0.0012 days, respectively; IRAS 18095+2704: +0.0027 and 
—0.0018 days, respectively; IRAS 19475+3119: +0.0021 and —0.0017 days, respectively. 
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Table 6 

Offsets Between the Published Data Sets and VUO Data Sets 11 


IRAS ID 

Source 

Dates 

Bandpass 

Offset (mag) 

17436+5003 

Fernie b 

1980-1998 

V 

0.000 


Percy & Welch c 

1979-1980 

V 

0.000 

18095+2704 

Arkhipova d 

1993-2008 

V 

-0.048 


Arklripova d 

1993-2008 

B 

-0.004 

19386+0155 

Arkhipova d 

1993-2008 

V 

0.060 


ASAS e 

2002-2009 

V 

-0.055 

19475+3119 

Arkhipova f 

2002-2005 

V 

-0.091 


a Offset = VUO — published data set. 
19861 


^ Fernie 

(19931 


(19831 


19891 1990al lbl. 1991 ): Fernie fe Seager 


19951) . 


Data from 1991 to 1998 are 


available at 
the other 


http: / /www. astro. utoronto. ca/DD 0 /research/89her .htrnl: 
data are listed in tables in the publications. 

c Percy fc Welch! ( 1981 ). where the data are listed. 

1 Arkhipova et al.l ( 19931 [2000l l2010t) : data available on-line from VizieR 
at http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=J/PAZh/26/705 and 
http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=J/PAZh/36/281. 

e littp://www.astrouw.edu.pl/asas/. 

‘ Arkhipova et al.l ( 20061 ): data available on-line from VizieR at 
http://vizier. u-strasbg.fr/viz-bin/VizieR?-source=J/PAZh/32/48 
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Table 7 

Periodogram Study of the Light Curves of Four O-Rich PPNs 1 


IRAS ID 

Set 

Filter 

Years 

No. 

Obs. 

Pi 

(days) 

Ai 

(mag) 

P 2 

(days) 

A2 

(mag) 

P 3 

(days) 

A 3 

(mag) 

P 4 

(days) 

a 4 

(mag) 

Ps 

(days) 

A 5 

(mag) 

Pe 

(days) 

A6 

(mag) 

17436 

A 

V 

1979-2012 

1640 

45.15 

0.014 

47.36 

0.015 

46.75 

0.012 

49.89 

0.014 

50.38 

0.012 

44.96 

0.013 

17436 

B 

V 

1979-1998 

1083 

45.08 

0.026 

46.54 

0.018 

47.45 

0.018 

49.58 

0.017 

55.47 

0.03 

39.62 

0.012 

17436 

C 

V 

1994-2012 

557 

47.30 

0.020 

49.92 

0.017 

42.36 

0.017 







17436 

D 

Rc 

1999-2012 

534 

47.36 

0.014 

49.82 

0.012 

42.34 

0.012 







17436 

E 

V 

1994-1999 

146 

44.09 

0.032 

42.38 

0.027 

48.90 

0.018 

72.43 

0.017 

53.89 

0.014 



17436 

F 

V 

2003-2007 

143 

41.47 

0.020 

56.87 

0.021 

43.88 

0.013 

54.15 

0.012 





17436 

G 

V 

2008-2012 

268 

47.42 

0.030 

57.70 

0.022 

38.14 

0.019 







18095 

A 

V 

1993-2012 

844 

113.3 

0.017 

97.9 

0.011 

103.0 

0.010 

108.4 

0.009 

158.9 

0.009 

116.1 

0.009 

18095 

B 

V 

1993-2002 b 

327 

98.2 

0.015 

107.5 

0.013 









18095 

C 

V 

2003-2012 

517 

111.8 

0.029 

110.2 

0.023 

103.1 

0.009 

68.3 

0.009 

97.2 

0.009 



18095 

D 

B 

1993-2012 

544 

113.5 

0.024 

103.0 

0.016 

154.5 

0.016 

97.9 

0.015 

116.3 

0.016 



18095 

E 

Rc 

1998-2012 

455 

102.9 

0.009 

112.6 

0.014 

109.1 

0.010 

97.7 

0.009 

149.3 

0.07 



18095 

F 

Rc 

2003-2012 

359 

101.1 

0.023 

99.9 

0.021 









19386 

A 

V 

1993-2008 

412 

101.8 

0.026 

98.6 

0.025 

103.5 

0.024 

106.3 

0.018 

89.2 

0.010 



19386 

B 

V 

1993-2000 

170 

99.5 

0.035 

91.7 

0.018 

109.3 

0.012 







19386 

C 

V 

2001-2008 

242 

103.2 

0.045 

98.8 

0.031 

106.7 

0.021 

88.3 

0.012 





19386 

D 

Rc 

2002-2007 

100 

100.4 

0.041 

105.1 

0.022 









19475 

A 

V 

1994-2012 

611 

40.95 

0.016 

38.85 

0.013 

48.09 

0.011 







19475 

B 

Rc 

1998-2012 

403 

40.98 

0.016 

38.84 

0.013 

48.18 

0.011 







19475 

C 

V 

2002-2012 

456 

40.99 

0.018 

38.79 

0.016 

48.11 

0.014 







19475 

D 

Rc 

2002-2012 

339 

41.02 

0.017 

38.80 

0.013 









19475 

E 

V 

1994-1999 

155 

36.72 

0.016 

36.13 

0.015 










a The formal, least-squares uncertainties determined using Period04 for each object in period P and amplitude A are approximately as follows: 

IRAS 17436: ±0.01 to ±0.025 days for subsets A to D and ±0.05 to ±0.20 days for the shorter subsets E to G, with ±0.001 to ±0.002 mag, respectively; 

IRAS 18095: ±0.1 to ±0.3 days for subsets A to E and ±0.8 days for subset F (strong beat periods), with ±0.001 to ±0.002 mag, except for subset F, in which 
they are much more uncertain; 

IRAS 19386: ±0.1 (subset A) to ±0.35 (subsets B to D) days, with ±0.002 mag; 

IRAS 19475: ±0.01 to ±0.04 days for subsets A to D and ±0.07 days for the shorter subset E, with ±0.002 mag, respectively. 
b A similarly good fit is achieved with Pi =98.5 days, Ai =0.013 mag and P2=102.8 days, A 2 = 0.012 mag 





















Table 8 

Results of the Period and Light Curve Study of Four Oxygen-Rich PPNs 


IRAS ID 

Pi 

(day) 

P2/P1 

p a 

^ prev 

(day) 

SpT 

Teff b 

(K) 

AV C 

(mag) 

A(V-R c ) b 

(mag) 

A(B—V) b 
(mag) 

AT efr b 

(K) 

Comments 

17436+5003 

45.2 

1.06 

40-46 

F3 lb 

6600, 7100 

0.21 

0.065 

0.095 

700 


18095+2704 

113.3 

0.86 

109 

F3 lb 

6500 

0.14 

0.060 

0.080 

590 

brightening trend 

19386+0155 

101.8 

0.97 

102, 98 

F5 lb 

6800 

0.20 

0.070 

0.105 

770 

brightening trend 

19475+3119 

41.0 

0.95 

43, 45 

F3 lb 

7750, 7200 

0.19 

0.045 

0.055 

400 



a Previousl y published perio ds: IRAS 17436+5003: F ernie, series of papers; IRAS 18095+2704: I Arkhipova et all (l2010h : IRAS 
19386+0155: I Arkhip ova et all (l2010h ; IRAS 19475+3119: lArkhipova et ahl (120061) . 

b T e ff determined from model atmosphere analyses as cited in Table [I] 

c The maximum brightness and color range observed in a season and the corresponding maximum temperature change based on 
A(B—V) and color—temperature calibration of lCoxI (|2000| , Table 15.7). 


to 

co 




















1995 2000 2005 2010 1994 1995 1996 199? 1998 1999 


1 17436+5003 -4V 04R-0.44 *4V(C1-C2)-3.B0 | 

-3.6 

M -3.5 

- R c „ k : 

:: ;• 1' 4* f : 

$ f l * 1 1 

£ $ i »ft}r *.<:< in'*?** 


-020 4 

+ ?• l| : 

50000 52000 54000 56000 

HJD - 2400000 


9500 50000 50500 51000 51500 

HJD - 2400000 

*• i* r % 1 «• 

“ - 3 - 5 

C? f & ft f 

i % 

-3.4 

-3.3 


* ^ | 9i fb 

-0.18 < 

% f : 


52000 52500 53000 53500 54000 54500 55000 55500 56000 

HJD - 2400000 HJD - 2400000 


Fig. 1.- Plots showing the differential light and 
color curves of IRAS 17436+5003 obtained at the 
VUO from 1993—2012. In the upper left panel are 
the entire curves and in the three other panels the 
curves are shown on expanded scales covering six 
or seven years each. The error bars are included 
with the data in these expanded plots. In the up¬ 
per left panel is also shown the differential V light 
curve of the primary (Ci) with respect to the sec¬ 
ondary (C 2 ) comparison stars on the same scale 
as the PPN light curves. The absence of V data 
from 2000 to part of 2002 is due to a filter prob¬ 
lem. Zero-point offsets are added to conveniently 
display all three light curves on the same plot. 
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Fig. 2.- Plots showing samples of the differen¬ 
tial B and (B— V) curves from the 2009 and 2010 
seasons obtained at the VUO. The changing cycli¬ 
cal behavior of the light and color curves can be 
seen, as can the correlation between brightness 
and color. Error bars are included. 



Fig. 3. Plots showing the change in color with 
change in brightness. A clear trend is seen, with 
the objects being redder when fainter. The ob¬ 
servations with the initial VUO CCD are plotted 
as filled circles and with the newer CCD as open 
circles. For IRAS 18095+2704, an adjustment in 
( V—Rc ) of 0.01 mag has been applied to correct 
for the small offset in color between the old and 
new CCD systems, as described in the text. 
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Fig. 4.— Plots showing the differential light and 
color curves of IRAS 18095+2704 obtained at the 
VUO, plotted similar to Figure [Q 


1994 1996 1998 2000 2002 2004 2006 2008 1994 1995 1996 1997 1998 


1 19386+0155 -iV "4R+0.10 4 4V(C1-C2)-1.42 | 

<■ -2.0 

- R ° * J* \ 

% 

>.*• w ,V 

1 * 

8 = «. V t $ $ 

. Al . j lii ;• = 

: v t i !? 

i ■:} /? A * 

'**•*?*#. ' * *• * ? *- 

* '' ■ * ‘!; f ’ i. t 

0i32 ‘f 
0.34 > 

i 

50000 51000 52000 53000 54000 

HJD - 2400000 


9500 50000 50500 51000 

HJD - 2400000 

. , v 


* ’ U.: it' 

- V * #5 i % 


Y * * ° f 

S 

-2.0 

!• f j! \i. 

V ; *V V 

> 


0.32 “ 

0-34 > 

•f # 1# 


51500 52000 52500 53000 53500 54000 54500 

HJD - 2400000 HJD - 2400000 


Fig. 5.— Plot showing the differential light and 
color curves of IRAS 19386+0155 obtained at the 
VUO from 1993—2007, plotted similar to Figured] 



' IT9475+3119 • 

V 04R+O.IO *4V(Cl-C2)-0.88 I " 





i*:> i- 1ft« <* 




•ij f;‘ V. K& 
a ’ # , 

* 

t. 

S. ■ 

B - " . 

, 

» 


• +■ 


/ * 

a*' 





’ ^ 

t 




Fig. 6.- Plot showing the differential light and 
color curves of IRAS 19475+3119 obtained at the 
VUO, plotted similar to Figure |H 
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Fig. 7. — Top panel: Combined Flight curve of 
IRAS 17436+5003, with our data shown as filled 
circles, the data of Fernie shown as open circles, 
and the data of Percy & Welch shown as filled 
triangles. Bottom five panels: The seasonally- 
adjusted V light curve fitted with the six peri¬ 
ods and amplitudes listed in Table [T] The fit is 
generally good throughout except for the most 
recent data (bottom panel). Note that there is 
some repetition of the data and the fitted curves 
at the edges of the five lower panels. The vertical 
dashed line shows the time when published high- 
resolution spectra were taken. 
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Fig. 8.- The frequency spectrum of the combined 
and seasonally-adjusted V light curves of the four 
targets, together with their associated phase plots 
based on the frequency peaks (Pi). The presence 
of secondary period peaks in the data is evident. 
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Fig. 9.— Top panel: Combined Flight curve of 
IRAS 18095+2704, with our data shown as filled 
circles and the data of Arkhipova et al. shown 
as open circles. Bottom two panels: The trend- 
removed, seasonally-adjusted V light curve on an 
expanded scale, fitted with the four periods and 
amplitudes listed in Table [7] The fit is reasonable 
on most seasons given the complex nature of the 
light curve. In some seasons, the periods looks 
like they agree but not the fixed amplitudes. The 
vertical dashed line shows the time when published 
high-resolution spectra were taken. 


1995 2000 2005 



Fig. 10. Top: Combined V light curve of IRAS 
19386+0155, with our data shown as filled circles 
and the data of Arkhipova et al. shown as open 
circles. Bottom: The trend-removed, seasonally- 
adjusted V light curve fitted with the five periods 
and amplitudes listed in Tabled] The fit is reason¬ 
ably good over this 16-year interval. The vertical 
dashed line shows the time when published high- 
resolution spectra were taken. 
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Fig. 11.— Top: Combined V light curve of IRAS 
19475+3119, with our data shown as filled circles 
and the data of Arkhipova et al. shown as open 
circles. Bottom three panels: The seasonally- 
adjusted V light curve fitted with the three pe¬ 
riods and amplitudes listed in Table [3 The fit is 
not so good, particularly in the amplitudes. We 
attribute this to the complex nature of the light 
curves, likely due to the presence of shock waves 
in the atmosphere. Note that there is some repeti¬ 
tion of the data at the edges of the two lower pan¬ 
els. The vertical dashed line shows the time when 
published high-resolution spectra were taken. 



Fig. 12.- A plot of pulsation period versus T e ff. 
The 12 C-rich objects a re shown with open cir¬ 
cles (Hrivnak et al.1 12010 ) and the four O-rich ob¬ 
jects from this study are shown with filled cir¬ 
cles. Those with less certain values are shown 
with smaller symbols. The straight line is a fit 
to the 12 C-rich objects. Two of the O-rich ob¬ 
jects (IRAS 17436+5003 and 19475+3119) have 
two points each because they have two different 
temperature measurements. (We have assumed 
reasonable error bars of ±250 K and ±3 days.) 
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Fig. 13.— Plots showing the comparison of the 
trends in the values of the maximum amplitude 
(AV) in a season with period (top) and tempera¬ 
ture (bottom). The open circles show the values 
for the 12 C-rich PP Ns and the dashed li nes show 
the general trends ( Hrivnak et al. 2010l) . Those 


with less certain values are shown with smaller 
symbols. The filled circles show the values for 
these four O-rich PPNs. The agreement is good 
between the two sets of data: they all show low 
amplitudes for shorter periods and higher temper¬ 
atures. (Error bars error bars of ±250 K, ±3 days, 
and ±0.02 mag.) 
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